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Q^' The Fortran code NMHDECAY computes the masses, couphngs and decay 

^ ■ widths of aU Higgs bosons of the NMSSM in terms of its parameters at the 

electroweak (or SuSY breaking) scale: the Yukawa couplings A and k, the soft 
trilinear terms Ax and A^, and tan/3 and fieS = A (S). The computation of 
^ ■ the spectrum includes leading two loop terms, electroweak corrections and 

propagator corrections. The computation of the decay widths is carried out 
as in HDECAY, but (for the moment) without three body decays. Each point 
in parameter space is checked against negative Higgs bosons searches at LEP, 
including unconventional channels relevant for the NMSSM. One version of 
the program uses generalized SLHA conventions for input and output. 
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1 Introduction 

The Next to Minimal Supersymmetric Standard Model (NMSSM [1-4,6-12]) 
provides a very elegant solution to the /i problem of the MSSM via the introduc- 
tion of a singlet superfield S. For the simplest possible scale invariant form of the 
superpotential, the scalar component of S acquires naturally a vacuum expectation 
value of the order of the SuSY breaking scale, giving rise to a value of fi of order 
the electroweak scale. The NMSSM is actually the simplest supersymmetric exten- 
sion of the standard model in which the electroweak scale originates from the SuSY 
breaking scale only. 

In addition, the NMSSM renders the "little fine tuning problem" of the MSSM, 
originating from the non-observation of a neutral CP-even Higgs boson at LEP II, 
less severe [2]. 

A possible cosmological domain wall problem [3] can be avoided by introducing 
suitable non-renormalizable operators [4] that do not generate dangerously large 
singlet tadpole diagrams [5]. 

Hence, the phenomenology of the NMSSM deserves to be studied at least as fully 
and precisely as that of the MSSM. Its particle content differs from the MSSM by the 
addition of one CP-even and one CP-odd state in the neutral Higgs sector (assuming 
CP conservation), and one additional neutralino. Thus, the physics of the Higgs 
bosons - masses, couplings and branching ratios [1,6-12] can differ significantly 
from the MSSM. The purpose of the Fortran code NMHDECAY (Non Minimal 
Higgs Decays), that accompanies the present paper, is an accurate computation of 
these properties of the Higgs bosons in the NMSSM in terms of the parameters in 
the Lagrangian. As its name suggests, the Fortran code uses to some extent - for 
MSSM-like processes - parts of the code HDECAY that is applicable to the Higgs 
sector of the MSSM [13]. 

In the present paper we define the NMSSM in terms of its parameters at the 
Susy breaking scale. No assumption on the soft terms (like universal soft terms at 
a GUT scale) are made. The parameters in the Higgs sector are chosen as follows: 

a) Apart from the usual quark and lepton Yukawa couplings, the scale invariant 
superpotential is 

XSHuH,+ '^S^ (1.1) 

depending on two dimensionless couplings A, k, beyond the MSSM. (Hatted 



capital letters denote superfields, and unhatted capital letters will denote their 
scalar components). 

b) The associated trilinear soft terms are 

XA^SH^H, + '^A^S' . (1.2) 

c) The final two input parameters are 

tan/5 = {Hu)l{Hd) , /Xeff = A (5) . (1.3) 

These, along with M^, can be viewed as determining the three SuSY breaking 
masses squared for iJ„, Hd and S through the three minimization equations 
of the scalar potential. 

Thus, as compared to two independent parameters in the Higgs sector of the MSSM 
(often chosen as tan/5 and M^), the Higgs sector of the NMSSM is described by the 
six parameters 

A , K , Aa , A^, tan/3 , /ieff . (1.4) 

We will choose sign conventions for the fields such that A and tan/3 are positive, 
while K, A\, Ak and fics should be allowed to have either sign. For any choice of these 
parameters - as well as of the values for the gaugino masses and of the soft terms 
related to the squarks and sleptons that contribute to the radiative corrections in the 
Higgs sector and to the Higgs decay widths - NMHDECAY performs the following 
tasks: 

1. It computes the masses and couplings of all physical states in the Higgs, 
chargino and neutralino sectors. Error messages are produced if a Higgs or 
squark mass squared is negative. 

2. It computes the branching ratios into two particle final states (including char- 
ginos and neutralinos — decays to squarks and sleptons will be implemented 
in a later release) of all Higgs particles. 

3. It checks whether the Higgs masses and couplings violate any bounds from 
negative Higgs searches at LEP, including many quite unconventional channels 
that are relevant for the NMSSM Higgs sector. It also checks the bound on 
the invisible Z width (possibly violated for light neutralinos). In addition. 



NMHDECAY checks the bounds on the hghtest chargino and on neutrahno 
pair production. Corresponding warnings are produced in case any of these 
phenomenological constraints are violated. 

4. It checks whether the running Yukawa couphngs encounter a Landau singu- 
larity below the GUT scale. A warning is produced if this happens. 

5. Finally, NMHDECAY checks whether the physical minimum (with all vevs 
non-zero) of the scalar potential is deeper than the local unphysical minima 
with vanishing {Hy) or {Hd)- If this is not the case, a warning is produced. 

The web sites 

http : //www . th . u-psud . f r/NMHDECAY/nmhdecay . html 
http : //higgs . ucdavis . edu/nmhdecay/nmhdecay . html 
provide links for downloading two different Fortran codes, along with descriptions, 
input files, and the data files that are required to apply the various LEP constraints. 

The Fortran code NMHDECAY_SLIIA.f reads from input files, and produces 
output files, that are as close as possible to the SuSY Les Houches Accord (SLHA) 
conventions [14]. Some generalizations of these conventions - including proposals for 
PDG numbers - have been necessary, however, in order to denote the NMSSM input 
parameters (1.4) and the additional particles in the Higgs and neutralinos sectors. 

The Fortran code NMHDECAY_SCAN.f reads from "private" input files (sam- 
ples are provided) that are constructed so as to scan over the NMSSM input parame- 
ters (1.4). Here the output is either "long" (easily human readable, if the number of 
points in parameter space is not too large), or "short", i.e. simple rows of numbers 
per point in parameter space, that should be edited according to the user's needs. 

Note that the "long" output also gives the reduced couplings of all neutral Higgs 
bosons to gauge bosons (CV), up type quarks (CU), down type quarks (CD), two 
gluons (CG) and two photons (CCA) (all relative to a Standard Model Higgs boson 
of the same mass). Using these, it is easy to compute the NMSSM Higgs production 
cross sections at colliders by rescaling those for the SM Higgs boson. 

The outline of the paper is as follows: In section 2, we describe in detail the 
accuracy with which the Higgs masses and mixing matrices are computed. In sec- 
tion 3, we describe the Higgs decay channels and corresponding accuracies that are 
used for the computation of the widths and branching ratios. In section 4, we de- 
scribe the various phenomenological constraints that can be applied to the model. 



Section 5 contains a user's guide for the Fortran codes, as well as our generalizations 
of the SLHA conventions. In section 6, we show some results for branching ratios 
produced with the present code, emphasizing decay channels that are atypical of the 
MSSM, such as the lightest CP-even state hi decaying to 2 lighter CP-odd states 
OiOi. Section 7 contains our conclusions. Appendix A contains a summary of our 
conventions, appendix B a summary of the Feynman rules, and appendix C the 
details of the radiative corrections to the higgs mass matrices. 

2 Radiative Corrections in the Higgs Sector 

Our convention for the superpotential, the soft terms, the resulting tree level 
potential and tree level mass matrices are given in appendix A. In the present 
section, we describe the accuracy with which radiative corrections are computed for 
the Higgs sector. 

First, we assume that the Yukawa couplings and soft terms are defined at a 
Susy breaking scale Q = Msusy, which corresponds to an average of the squark 
masses. Quantum fluctuations at scales > Q are assumed to be integrated out, which 
corresponds to the standard RG evolution of the parameters from a fundamental 
scale like Mqut down to Msusy- The effective Lagrangian at the scale Q can be 
assumed to be of the standard supersymmetric form plus soft terms. 

We are interested in the full effective action 

Feff = J2Z, D^H^D^H, - \/eff(//.) , (2.1) 

i 

that is obtained from the effective Lagrangian at the scale Q by adding all quantum 
effects arising from fluctuations at scales < Q. (Here, Hi denotes all the Higgs flelds, 
Hu, Hfi and S.) These quantum effects can be classifled according to powers of the 
various couplings, and powers of potentially large logarithms. Explicit formulas 
for the radiative corrections to Vg//, Zi and hence to the Higgs mass matrices are 
given in appendix C, subsequently we describe only the orders to which radiative 
corrections are taken into account. 

Let us start with the corrections to the effective potential. It is somewhat more 
convenient, however, to classify the corrections to the (lightest) scalar Higgs mass, 
which is essentially the second derivative of the effective potential. At tree level this 



mass squared is bounded by 



m\<Ml\ cos^ 2/^ + — 2 2^ ) , (2.2) 



which is of the order 



m\ ~ C((?2t;2) + Oiy^v") (2.3) 

where (7 denotes the electroweak gauge couphngs, and v the magnitude of the Higgs 
vacuum expectation values (we do not distinguish between large and small tan /3 in 
this section). 

The dominant one-loop corrections to vn\ originate from top, stop, bottom and 
sbottom loops. The corresponding corrections bm\ to vn\ are of order 

f^top"^? ~ /^t v^ In (QV"^?op) , 

4top"^? ~ h\ v^ (for g2 _ ^2^^p) ^ 

(^bottomm^ ~ /l^ v^ In (QV^^bottom) > 

4bottom"^? ~ hi V^ (for Q2 ~ "^sbottom) • (2-4) 

We have included these contributions [6,7, 11] exactly, without expanding in large 
logarithms or squark mass splittings. 

The dominant two loop corrections are of the form [11, 12] 

5(2)m? ~ hy In^ (gV^p) + h\a,v^ In^ (QV^p) ■ (2-5) 

Here, we only include the leading double logarithms, i.e. we neglect single logs as 
well as terms involving /i^. 

Coming back to the one-loop order, we focus on including all contributions of the 
order g^^ large logarithms. These can be obtained by RG methods [15], taking care 
of the fact that the electroweak couplings g are defined at the scale Mz and not at the 
scale MsusY- We have included all these contributions as given in [15], distinguishing 
large logarithms ln(g7M|) from ln(M|/M|) and ln(/x7M|) (/i^ = Mi,M2 or /ieg, 
where Mj are the electroweak gaugino masses). We have not (yet) included terms 
of the order g'^X^x large logarithms, nor terms of the order (yf^x finite. 

This concludes our contributions to the effective potential. Next, we outline the 
contributions to the wave function normalization constants Zi in (2.1). If evaluated 
for external momenta of (9(mtop) (the order of the Higgs masses), top and bottom 
quark loops yield contributions to 5Zi (relative to 1) of the form 

5Z, ~ hi In (QV^Lp) , 5Z, ~ hi In (QV^Lp) • (2-6) 



After rescaling the Higgs fields so that their kinetic energies are canonically normal- 
ized, these effects generate contributions dzml which take the form 

5zml r. gX/,v^ In (gV^Lp) + \X/y In (QV^p) • (2-7) 

However, after diagonalization of the Higgs mass matrices, one may find eigen- 
states with masses rriH substantially larger than mtop- Then the Z factors should 
be evaluated at external momenta of 0{mH), i-e. at the pole of the propagators. 
Hence for rriH ^ mtop the logarithms lia^Q"^ / ml^p) in (2.6) should be replaced by 
\n(Q'^/m'jj), with coefficients depending on Higgs mixing matrices. 

All this is done in the program, with the net result that all potentially large 
logarithms with coefficients g'^h'^i^, or A^/iwj^ have the correct arguments. (We neglect, 
however, contributions of 0{\n{m'\j/m^ )) in the off diagonal Higgs matrix elements, 
which affects only the couplings of very heavy Higgses.) Since we have neglected 
terms of 0{h^n x finite) in the Z factors and electroweak corrections without large 
logarithms, the dominant sources of uncertainty in Sm^ are contributions of 0{g'^v'^ x 
finite), 0{g'^hfnv'^ x finite) and 0{\^h1nv'^ x finite) as well as subdominant (single 
log) two loop corrections. 

We have checked that in the end, after expanding the Z factors appropriately, 
all Q^ dependence in the one-loop contributions can be reabsorbed into the running 
of the input parameters that are defined at the scale Q^. An exception is the Q^ 
dependence of terms ~ g^, since we have not included the squark loop contributions 
~ g^ (from squark couplings to Higgs bosons arising from the D-terms). These 
would give contributions to m\ of the order 

5Dml^gHn{Q'/ml^^.^,^) , (2.8) 

which is relatively small for Q^ ~ ""^squark- 

Since we have not included contributions of the orders {g^, g'^hfi^)x finite, we are 
actually not sensitive to the scheme in which the electroweak couplings are defined. 
Since in the computation of the decay widths we use the physical masses Myy and 
Mz in order to parameterize some couplings, the most useful scheme here is the on- 
shell scheme, where the three "unphysical" electroweak parameters f ^ = {H^ + -ffj). 



gl and (/I are defined in terms of Gp, Mw and Mz as 

1 



v^ 



2V2Gf ' 
gl = 2M^/v\ 

gl + gl = 2M|y. (2.9) 

Finally we remark that we have applied the above corrections consistently to the 
complete Higgs sector, i.e. the CP-even, CP-odd and charged Higgs states. One finds 
that one can reabsorb several of the radiative corrections into a redefinition of the 
input parameter A\ [6], after which the CP-odd mass matrix assumes its tree level 
form, up to the required rescalings of the vevs by the Z factors. This phenomenon 
is familiar from the MSSM, where Ma is thus a convenient input parameter, and 
the required rescalings are often absorbed into a scale dependent value of tan (3. 

We have checked that in the MSSM hmit of the NMSSM (A, /t < 1), after 
comparing for the same values of tan/? and Ma (and the remaining soft terms), the 
mass of the lightest CP-even state agrees with the one computed in HDECAY [13] 
(models 1 [16] or 2 [18]) to better than 2 GeV for moderate mixing (A^p < 1 TcV), 
and to better than 3 GeV for maximal mixing. This coincides with the expected 
theoretical error on this mass given the orders of uncertainty discussed above. 

A more detailed discussion of the radiative corrections to the NMSSM Higgs 
sector will appear in [17]. 

3 Higgs Decays 

In this section we describe the decay modes that are included in the Fortran 
code. Most of the corresponding code is extracted from HDECAY [13]. In what 
follows, H denotes any of the 3 CP-even or 2 CP-odd scalars, or the charged Higgs 
boson. 

a) H ^ gluons: We take into account charm, bottom and top quark loops; 
the lowest order contribution to the decay width is given in [19]. QCD radiative 
corrections [20] are included at the same level as in HDECAY. Squark loops are not 
(yet) taken into account. 

b) iJ — * leptons and quarks: QCD corrections are included as in HDECAY [21], 
but we use always the two loop formula for the relation between the pole quark mass 
mpoic and the running MS quark mass m{mpo\c)'- 



mi{mpo\e) = ^^^^— -2 (3.1) 

1 _l_ 4as(mpoic) ry- / as(mpoic) \ 

Stt V ■"" / 

with ^bottom ~ 12.4, Ktop ~ 11- The scale dependence of the running s and c 
quark masses is computed as in HDECAY, but for the running b and t quark masses 
we use running Yukawa couphngs that solve the one-loop RG equations exactly in 
as, improved by corrections to first order in the Yukawa couplings themselves (see 
appendix C). 

Since electroweak corrections to decays into leptons and quarks are small [13], 
they are neglected. Below-threshold three-body decays H ^ tt ^ tbW~ are not 
(yet) taken into account. However, contributions to the decay rates into hb {cc) 
from H -^ gg* -^ gbb {cc) are included as in HDECAY (option NF-GG = 3, that 
can easily be modified by re-editing the line NFGG = 3 at the beginning of the 
subroutine DECAY). SuSY loop corrections are not (yet) included. (These could be 
important for light squarks and/or large tan/?.) In the case of a decaying charged 
Higgs boson, we must input the CKM matrix elements V^s, V^b and Vcb- 

c) H ^ WW/ZZ: We treat these decays as in HDECAY [22], but without the 
double off-shell option. 

d) H ^ 77/Z7: We take into account loops of charm, bottom and top quarks, r 
leptons, W bosons, charginos and charged Higgs bosons [23]. Loops of squarks and 
sleptons are not (yet) considered. 

e) H ^ HH: All kinematically possible combinations {H = CP-even, CP-odd 
or charged) are considered, with (Higgs)^-couplings computed from the NMSSM 
potential (see F. Franke and H. Fraas in ref. [1] and appendix B). The dominant 
(leading log) radiative corrections from top-quark and bottom-quark loops to these 
couplings are included in the Fortran code, the corresponding formulas will appear 
in [17]. Decays into three particle final states are not (yet) taken into account. 

i) H ^ HZ/ HW: All kinematically possible combinations {H = CP-even, CP- 
odd or charged) are considered [24]. Below threshold decays into three particle final 
states are not (yet) taken into account. 

g) H -^ charginos/neutralinos: First, the 2x2 chargino and 5x5 neutralino mass 
matrices are diagonalized in the subroutines CHARGINO and NEUTRALINO, and 
the decays into all possible two body final states [25] are included in the program. 

Decays into sleptons or squarks are not (yet) considered. Note that relatively 



light squarks - even if too heavy to be produced in Higgs decays - could nevertheless 
affect notably the processes b) and d) above through squark loops. For this reason 
the present version of NMHDECAY should not be applied to scenarios with light 
squarks. 

4 Constraints on the Parameter Space 

Once the spectrum and decay widths of the Higgses are computed, we check 
some experimental and theoretical constraints for each point in parameter space: 

First, we require that there is no Landau pole for the Yukawa couplings A, k,, /itop 
and /ibot below the GUT scale. For this, we integrate the two loop renormalization 
group equations for the Yukawa and gauge couplings from the SuSY scale up to the 
point where the gauge couplings gi and g2 unify. At this point we check that none of 
the couplings is too large (i.e. exceeds v^Itt). (For simplicity we assume a universal 
Susy threshold at Msusy ~ Q-) 

We also check that there is no deeper minimum of the scalar potential with 
vanishing (i/„) or {H^). 

Finally we check all available experimental constraints from LEP: 

1) In the neutralino sector, we check that the lightest neutralino does not contribute 
too much to the invisible Z width {T{Z -^ x?X?) < 1-76 MeV [26]) if m^o < Mz/2, 
and that o-(e+e~ -^ XiX^) < 10"^ pb iim^o +m^o < 209 GeV {i > 1) and o-(e+e~ -^ 
Xi'Xj) < 10^^ pb if fn^a + m^o < 209 GeV {i,j > 1) [27]. In the chargino sector, we 
check that the lightest chargino is not too light {iTi + < 103.5 GeV [28]). 

2) In the charged Higgs sector, we check the bound 'mH+ > 78.6 GeV [29]. 

3) In the neutral Higgs sector, we check the constraints on the production rates 
(reduced couplings) x branching ratios versus the masses, for all of the CP-even 
states h and CP-odd states a, in the following channels studied at LEP: 

• e"'"e^ -^ hZ, with h ^ bb and h -^ t^t^ (from the LEP Higgs working group 
results [30]); 

• e^e~ -^ hZ, with h decaying to two jets. For this, we combined the low mass 
range results from OPAL at LEP2 [31] with the higher mass range obtained 
at LEP2 by the LEP Higgs Working group 



e^e -^ hZ ^ with /i -^ 77 (from the LEP Higgs working group results [33]); 



10 



• e^e' -^ hZ^ with h decaying invisibly {i.e. into two neutralinos) . For this, 
we combined the low mass range results from ALEPH at LEPl [34] with the 
higher mass range obtained at LEP2 by the LEP Higgs Working group [35]. 

• e^e' -^ hZ, independent of the h decay mode, looking for a peak of the Mx 
recoil mass distribution in e+e^ -^ XZ. For this, we combined the low mass 
range results from ALEPH at LEPl [34] with the higher mass range obtained 
at LEP2 by OPAL [36]. 

• For the associated production mode e^e^ -^ ha we used the DELPHI results 
[37] for the various final states: ha -^ 46, ha -^ At and ha -^ aaa -^ 66. 
The latter decay channel plays an important role in the NMSSM, where the 
lightest CP-even Higgs can decay mainly into two light CP-odd Higgses over 
large areas of the parameter space, as we will see in section 6. (In the context 
of the CP-conserving MSSM, this is possible only for very special parameter 
choices.) DELPHI also studied the channel hZ -^ aaZ -^ 4b + 2jets [37]; 
their limits have also been implemented in our phenomenological constraints. 

• The channels e~^e^ -^ hZ -^ aaZ with aa -^ 4jets, 2jets + cc, 2jets + tt, 
4r's, cccc, TT + cc have been studied by OPAL [38] (see their Fig. 7 — we 
employ numerical files for the contours provided by the authors). These limits 
have also been included. 

All these constraints are contained in the files located in the directory LEPCON 
downloadable together with the Fortran codes. The file names are constructed so 
that they can be identified with the above delineated constraints. (The current 
numerical limits contained in these files can possibly be improved with the help of 
further dedicated analysis of existing data.) 

5 How to use NMHDECAY 

NMHDECAY exists in two versions: 

1. NMHDECAY_SLHA uses an input file and produces output files that are suit- 
able generalizations of the SLHA conventions [14]. It is designed for studying 
the properties of one user-defined point in parameter space. 
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2. NMHDECAY_SCAN uses input and output files that can be modified by the 
user so as to scan over parts of, or all of the NMSSM parameters defined in 
eq. (1.4). 

Both programs have the common feature that they are based on one single Fortran 
code (NMHDECAY_SLHA.f or NMHDECAY_SCAN.f) that does not need to be 
linked with any other code. However, both programs need data files in order to 
check against negative Higgs searches at LEP in the numerous channels discussed 
in section 4. These data files are available in the directory LEPCON, that can be 
downloaded from the NMHDECAY home page. Note that the directory LEPCON 
has to be situated in the same directory that contains the executable NMHDECAY 
code. We now discuss the particular features of the two programs. 

5.1 NMHDECAY_SLHA 

NMHDECAY_SLHA uses an input file slhainp.dat, a version of which is down- 
loaded automatically with the Fortran code. This sample file appears in Table 1. 
Several comments on its contents are in order. 

a) "BLOCK MODSEL" contains the entry 3 (corresponding to the choice of the 
particle content) with switch 1, as attributed to the NMSSM in [14]. 

b) "BLOCK SMINPUTS" contains important Standard Model parameters. 

1. First, there is the inverse electromagnetic coupling constant at the scale 0, 
which is that required for the computations of the decay widths into two (on- 
shell) photons. 

2. Second, various Higgs couplings are defined in terms of Mz-, as{Mz), My/ and 
Gpi which means an on shell scheme is implicitly being used in order to define 
the electroweak parameters (cf. the corresponding discussion in section 2). 

Mz, as{Mz) and Gp are input as part of "BLOCK SMINPUTS", whereas 
the numerical value of Myy = 80.42 GeV is defined in the subroutine INPUT, 
that reads also the input file slhainp.dat. (Part of the corresponding code 
is copied from pyslha.f , with thanks to P. Skands.) 

3. Third, as part of this block we input the running b quark mass mi,{mh), the 
top quark pole mass and m,-. 
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4. In addition, NMHDECAY needs the (pole) quark masses rris and rric, as well 
as the CKM matrix elements Vus, Vcb and Vub- The numerical values of these 
five parameters are defined in the subroutine INPUT. (For convenience, they 
are printed out in the output file spectr.dat, see below.) 

c) "BLOCK MINPAR" contains only the switch 3 with the input value for tan/3 
specified. 

d) "BLOCK EXTPAR" contains inputs for the SuSY and soft-SuSY-breaking pa- 
rameters. Needed generalizations of the SLHA conventions appear. The new entries 
are: 



61 


for A 


62 


for K 


63 


for Aa 


64 


for A,, 



Note that the entry 23 {ji) is used for the effective /i parameter, ii^a = A (S'), in the 
NMSSM. It should also be noted that neither the slepton entries (all generations) 
nor the squark entries for the first and second generations are actually used; these 
are included for future use only. 

The two output files of NMHDECAY_SLHA are spectr.dat (see Table 2) and 
decay.dat. In spectr.dat, "BLOCK SPINFO" is followed by warnings (switch 
3) if any of the constraints described in section 4 are violated, or if any of the six 
Higgs states could decay into squarks in which case the corresponding widths and 
branching ratios should not be used. This segment of the output also displays error 
messages (switch 4) if any of the Higgs, stop or sbottom states have a negative mass 
squared. No spectrum output is produced in this case. 

In the output, "BLOCK SMINPUTS" is followed by a printout of the Standard 
Model input parameters. The numerical values for Mw-, iTig, rric, Vus, Vcb and Vub, 
that have no SLHA numbers, appear in lines subsequent to "7^ SMINPUTS Beyond 
SLHA" . Similarly, "BLOCK MINPAR" is followed by a printout of the value of 
tan/9 and "BLOCK EXTPAR" is followed by a printout of the important SuSY and 
soft-SuSY-breaking parameters. 

The (pole) masses for the Higgs particles, the neutralinos, the charginos and 
squarks appear in the output following "BLOCK MASS". There, one finds several 
essential NMSSM generalizations of the SLHA conventions. The new entries, with 
proposed PDG codes, are 
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45 for the third CP-even Higgs boson, 

46 for the second CP-odd Higgs boson, 
1000045 for the fifth neutrahno. 

The Higgs mixings in the CP-even sector follow "BLOCK HIGMIX" and those in 
the CP-odd sector follow "BLOCK AMIX". Both segments of the output contain 
NMSSM generalizations of the SLHA conventions, that are required in order to 
parameterize the mixing in the enlarged Higgs sector. The meaning of the matrix 
elements Sij (i, j = 1, 2, 3) and P[j (i, j = 1, 2) is as follows. 

• According to the SLHA conventions the Higgs states Hu-, H^ are denoted by 
H2-, Hi, respectively. (Inside the Fortran code the Higgs states Hu-, H^ are 
denoted by Hi, H2, which is of no relevance for the SLHA output.) Hence, 
for the purpose of the SLHA output, the CP-even Higgs weak eigenstates are 
numbered by S'f ^"'^ = {HdR, H^r, Sr) {R refers to the real component of the 
field). If hi are the mass eigenstates (ordered in mass), the convention is 

I, c Qweak 

Ibi — OijOj 

• In the CP-odd sector the bare Higgs fields are H^i, H^i, Sj (/ for imaginary 
component). Again, for the purpose of the SLHA output, the CP-odd Higgs 
weak eigenstates are denoted by Hui = H21, H^i = Hu. The mass eigenstates 
are Oj (ordered in mass, i = 1,2) and the Goldstone mode G. Then the 
elements of P/,- are defined as 

Hu = smP{P[iai + -P2i'^2) + cos/^G , 
H21 = cosP{P[iai + ^21*^2) — sinPG , 
Si = P{2ai + P22a2 ■ (5.1) 

• In the output, "BLOCK NMIX" is followed by a printout of the obvious gen- 
eralization of the 4x4 MSSM neutralino mixing matrix to the 5x5 NMSSM 
neutralino mixing matrix (with real entries); "BLOCK UMIX" and "BLOCK 
VMIX" are followed by printouts of the U and V matrices as defined in the 
MSSM. 

• The output file decay.dat respects the SLHA conventions for the Decay file 
(at present we consider two particle final states only), using the above gen- 
eralizations of the PDG codes both for the decaying particle and the final 

states. 
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5.2 NMHDECAY_SCAN 

NMHDECAY_SCAN uses the input files smpar.dat and scaninp.dat, versions of 
wliicli are downloaded automatically with the Fortran code. One of the sample files 
appears in Table 3. 

In the file smpar . dat the following standard model parameters must be specified: 
as{Mz), Gf, c^iTm.lO)) ^^e lepton masses rrir and m^, Mz, Mw, the pole quark masses 
TTis, rric, the running bottom quark mass mi){mb), the top quark pole mass rrit, and 
the CKM matrix elements Vus, Vcb and Vub- 

In the file scaninp.dat, the following must be specified: 

• the total number of points to be scanned in parameter space; 

• the output format (0 for "short" , corresponding to simple rows of numbers per 
allowed point in parameter space, and 1 for "long", as described below); 

• lower and upper limits for the NMSSM parameters A, /t, tan/3, /ieg, A\ and 
^k; 

• the soft squark masses, trilinear couplings and gaugino masses over which no 
scan is performed — the slepton masses as well as the squark masses of the 
first two generations (without index 3) are not (yet) used. 

The scan in parameter space uses a random number generator, such that all NMSSM 
parameters are randomly chosen point by point in the parameter space within the 
specified limits. 

The output file containing the physical parameters is always called scanout . dat, 
regardless of the output format chosen. The numbers printed out for the output 
format (recommended for scans over more than 10 points in parameter space) 
should be edited according to the user's needs. ^ The output format 1 is easily 
readable and shows 

- the NMSSM parameters for each point, 

- possible warnings in case any phenomenological constraint is violated, a strong 
coupling regime below Mqut appears or an unphysical minimum is deeper than the 
physical one (cf. section 4), if any of the six Higgs states could decay into squarks 
(in which case the corresponding widths and branching ratios should not be used) 



^See the section of the program following the comment line " The following 3 lines should be 
edited according to the user's needs" . 
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or error messages ("fatal" errors) in case any of the Higgs, stop or sbottom states 
has a negative mass squared (in which case no additional output is produced). 

- for each of the six Higgs states, their mass, their decomposition into weak eigen- 
states {Hu, Hd, S) (in Components), their reduced couplings to gauge bosons (CV), 
up type quarks (CU), down type quarks (CD), two gluons (CG) and two photons 
(CGA) (all relative to a standard model Higgs boson with the same mass), their 
branching ratios (where "Higgses" denote all possible two Higgs final states, and 
"sparticles" all possible two particle neutralino/chargino final states), and their to- 
tal width, 

- the neutralino and chargino masses (all masses in GeV), as well as the neutralino 
composition in the basis ip^ = {—i\i,—i\2,ipuy'^dy'^s) described in detail in Ap- 
pendix A. 

- the (pole) stop and sbottom masses (in GeV). 

The output file scanerr.dat shows how many of the points in parameter space 
have avoided fatal errors or violations of phenomenological constraints, and the 
range in the NMSSM parameter space over which points have passed all these tests. 

For users who wish to call a subroutine as a function of the Higgs, chargino, 
neutralino and squark sector outputs, including mixing angles and other parameters 
and quantities computed during the course of the scan, they should use the param- 
eters and common blocks found in NMHDECAY_SCAN.f within "SUBROUTINE 
OUTPUT" . The comments included in this subroutine should allow easy identifica- 
tion of all the parameters, branching ratios, mixing angles and so forth that would 
be of potential interest for inputting into a user's subroutine. 

6 Results and discussion 

As stated above, the masses, couplings and decay properties of the Higgs bosons 
of the NMSSM can differ significantly from the MSSM. The primary purpose of 
NMHDECAY is to allow for detailed studies of such cases. 

A particularly interesting case is where the lightest Higgs scalar hi, although 
primarily non-singlet, decays mainly into two pseudoscalars (also primarily non- 
singlet) [9]. It is then interesting to see how the branching ratios of hi depend on 
the parameters of the model. In particular, we can determine whether the choices 
of parameters that yield the above types of decays are ruled out for some of the 
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phenomenological reasons discussed in section 4. We can also determine if there are 
particular (fine-tuned) relations between the parameters required for hi -^ aiai to 
be dominant. To study such issues, it will be convenient to fix all but one of the 
parameters (which allows for a reasonable graphical representation), and perform a 
scan over the remaining parameter. 

In some sense, the input parameter Ax is the most natural one to vary, since the 
mass of the MSSM like pseudoscalar depends quite strongly on Ax (and hence Ax 
plays the role of Ma in the MSSM). 

Let us first consider the following choice of the NMSSM parameters [cf. eq. (1.4)]: 
X = K, = 0.3, tan/? = 5, fi^s = 180 GeV, A^ = 0. For the squark masses and trilinear 
couplings, we take 1 TeV and 1.5 TeV, respectively. Varying A\ between and 
1000 GeV, we obtain the branching ratios for hi as shown in fig. 1. These show 
clearly that, for Ax < 600 GeV, the decay hi -^ ai Oi is dominant. 

The reason for the sharp drop of the hi -^ aiai branching ratio for Ax > 
600 GeV becomes clear from fig. 2, where we show the masses rrih^ and rria^ as 
functions of A\; for Ax > 600 GeV, rriaj^ becomes larger than m/i^/2. The hi -^ aiai 
decay is also reduced as A\ -^ 0, even though the ai becomes very light, because in 
this case it is mainly singlet and has a very small coupling to the hi. None of the 
points in these two graphs are excluded by LEP. They should be visible at the LHC 
using the techniques we have developed [9] for isolating the WW ^> h ^> aa type of 
signal. 

For a second sample set of plots, figs. 3-6, we take A = 0.5, k = —0.15, 
tan/5 = 3.5, fi^s = 200 GeV, Ax = 780 GeV and A^ e [150 GeV, 250 GeV]. The 
scaninp . dat file for this case is given in Table 3. For much of this parameter range, 
neither the hi nor the /i2 would have been observable at LEP. In particular, fig. 5 
shows that mfi2 ^120 GeV implying that the /i2 is beyond the LEP kinematical 
reach. The hi is much lighter. However, this light Higgs is not excluded by LEP 
over most of the above A,^ range since: a) its reduced coupling to gauge bosons 
is small; and b) hi -^ bb is suppressed so that hi -^ jj decays are dominant (see 
fig. 3). In fig. 6, we plot ^^ = Cv(hiY x BR{hi -^ jj) for our selected points as 
well as the region excluded by LEP searches in this channel [32]. We see that only 
if "^?ii ^ 53 GeV, which corresponds to A^ > 235 GeV, would the hi be excluded 
by LEP data. 

Will these Higgs bosons be observable at the LHC? In this regard, it is important 
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to note from fig. 4 that when A^ > 215 GeV, /i2 -^ hihi decays are dominant. This 
occurs because rrih^ decreases with A^, see fig. 5. Meanwhile, fig. 3 shows that 
BR{hi -^ bb) and BR{hi -^ t+t-) are both small when A^ e [205 GeV, 220 GeV]; 
in this region of parameter space, the hi decays mainly to cc or gg. Thus, for 
A^ ~ 215 - 220 GeV: 

• The hi has a mass that lies below the mass range currently studied for Higgs 
detection at the LHC. Further, the hi will be so weakly produced at the LIfC 
(since C,"^ < 0.1) that extensions to lower Higgs masses of the current LHC 
studies would probably conclude it was undetectable. 

• Simultaneously, the strongly produced /i2 has decays dominated by /i2 — * hihi 
with hi -^ cc,gg (but not bb or r^r^). As a result, the techniques of [9] for 
h ^ aa (which require a significant a -^ r^r^ branching ratio) do not apply, 
and the /i2 would also appear to be very difficult to observe at the LHC. 

We note that the above choice of parameters produces a phenomenology for the 
Higgses that is somewhat similar to that discussed in [10], which focuses on a region 
of parameter space such that the hi has suppressed decays to bb and r^r^ . The 
difference is that the decay of their /i2 to hihi is kinematically forbidden and the /i2 
could be detected at the LHC. 

7 Conclusions 

In this paper, we have presented details regarding the now publicly available 
NMHDECAY programs, NMHDECAY_SLHA.f and NMHDECAY_SCAN.f, that can 
be used to explore the Higgs sector of the Next to Minimal Supersymmetric Model 
defined by adding one singlet Higgs superfield to the Minimal Supersymmetric 
Model. The programs, and associated data files, can be downloaded from the two 
web pages: 

http : //www . th . u-psud . f r/NMHDECAY/nmhdecay . html 
http : //higgs . ucdavis . edu/nmhdecay/nmhdecay . html 
The web pages provide simplified descriptions of the programs and instructions on 
how to use them. The programs will be updated to include additional features 
and refinements in subsequent versions. We welcome comments with regard to 
improvements that users would find helpful. 
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We have exemplified the use of these programs for two particularly interesting 
scenarios: The first illustrates the potentially crucial importance of the LHC h -^ 
aa detection mode (that is dominant over a significant, not fine-tuned, range of 
parameters of the NMSSM). The second exposes a limited portion of parameter 
space for which Higgs discovery would not have occurred at LEP and will probably 
not be possible at the LHC. 
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# INPUT FILE FOR MMHDECAY 

# BASED ON SUSY LES HDUCHES ACCORD, MODIFIED FOR THE NMSSM 

# IN EXTPAR: LINES 61-64: NMSSM YUKAWA COUPLINGS AND TRILIN. SOFT TERMS 
BLOCK MDDSEL 

3 1 # NMSSM PARTICLE CONTENT 
BLOCK SMINPUTS 



1 


137.036 


# 


ALPHA_EM--1(0) 


2 


1.16639D-! 


5 # GF 


3 


0.12 


# 


ALPHA_S(MZ) 


4 


91.187 


# 


MZ 


5 


4.24 


# 


MB (MB), RUNNING B QUARK MASS 


6 


175. 


# 


TOP QUARK POLE MASS 


7 


1.7771 


# 


MTAU 


BLOCK 


MINPAR 






3 


5. 


# 


TANBETA 


BLOCK 


EXTPAR 






1 


5.D2 


# 


Ml 


2 


1.D3 


# 


M2 


3 


3.D3 


# 


M3 


11 


1.5D3 


# 


ATOP 


12 


1.5D3 


# 


ABOT 


13 


1.5D3 


# 


ATAU 


23 


180. 


# 


MU 


31 


1.D3 


# 


LEFT SELECTRON 


32 


1.D3 


# 


LEFT SMUON 


33 


1.D3 


# 


LEFT STAU 


34 


1.D3 


# 


RIGHT SELECTRON 


35 


1.D3 


# 


RIGHT SMUON 


36 


1.D3 


# 


RIGHT STAU 


41 


1.D3 


# 


LEFT 1ST GEN. SQUARKS 


42 


1.D3 


# 


LEFT 2ND GEN. SQUARKS 


43 


1.D3 


# 


LEFT 3RD GEN. SQUARKS 


44 


1.D3 


# 


RIGHT U-SQUARKS 


45 


1.D3 


# 


RIGHT C-SQUARKS 


46 


1.D3 


# 


RIGHT T-SQUARKS 


47 


1.D3 


# 


RIGHT D-SQUARKS 


48 


1.D3 


# 


RIGHT S-SQUARKS 


49 


1.D3 


# 


RIGHT B-SQUARKS 


61 


.3D0 


# 


LAMBDA 


62 


.3D0 


# 


KAPPA 


63 


200. 


# 


A_LAMBDA 


64 


0.0 


# 


A_KAPPA 



Table 1: Sample slhainp.dat file. 
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# NMHDECAY OUTPUT IN SLHA FORMAT 

# Info about spectrum calculator 



BLOCK SPINFO 

1 NMHDECAY 

2 1.1 

# Input parameters 
BLOCK MODSEL 

3 1 
BLOCK SMINPUTS 

1 

2 

3 



# Program information 

# spectrum calculator 

# version number 



# NMSSM PARTICLE CONTENT 



1.37036000E+02 
1.16639000E-05 
1.20000000E-01 

4 9.11870000E+01 

5 4.24000000E+00 

6 1.75000000E+02 

7 1.77710000E+00 
# SMINPUTS Beyond SLHA: 



# ALPHA_EM~-1 

# GF 

# ALPHA_S(MZ) 

# MZ 

# MB (MB) 

# MTOP (POLE MASS) 

# MTAU 



# MW 

# MS 

# MC 

# VUS 

# VCB 

# VUB 



0.80419998E+02 
0.19000000E+00 
0.16100000E+01 
0.22200000E+00 
0.41000000E-01 
0.36000000E-02 



BLOCK MINPAR 
3 5. 
BLOCK EXTPAR 



OOOOOOOOE+00 # TANBETA 



1 
2 
3 
11 
12 
13 
23 
33 
36 
43 
46 
49 
61 
62 
63 
64 



5.00000000E+02 
l.OOOOOOOOE+03 
3.00000000E+03 
1.50000000E+03 
1.50000000E+03 
1.50000000E+03 
1.80000000E+02 
l.OOOOOOOOE+03 
l.OOOOOOOOE+03 
l.OOOOOOOOE+03 
l.OOOOOOOOE+03 
l.OOOOOOOOE+03 
3.00000000E-01 
3.00000000E-01 
2.00000000E+02 
0. OOOOOOOOE+00 



# Ml 

# M2 

# M3 

# ATOP 

# ABOTTOM 

# ATAU 

# MU 

# LEFT STAU 

# RIGHT STAU 

# LEFT 3RD GEN. S QUARKS 

# RIGHT T-SQUARKS 

# RIGHT B-SQUARKS 

# LAMBDA 

# KAPPA 

# A_LAMBDA 

# A KAPPA 



# 

BLOCK MASS 

# PDG Code 

25 

35 

45 

36 

46 

37 

1000022 

1000023 

1000025 

1000035 

1000045 

1000024 

1000037 

1000005 

1000006 

2000005 

2000006 



# Mass spectrum 
mass 

1.16801961E+02 
3.56992533E+02 
5.91271441E+02 
4.90592823E+01 
5.87736637E+02 
5.90714771E+02 
1.66456703E+02 
-1.85971845E+02 
3.67762262E+02 
5.04634098E+02 
1.00711878E+03 
1.76259627E+02 
1.00710838E+03 
8.21132887E+02 
6.24445725E+02 
8.24131926E+02 
1.00488215E+03 



particle 

# lightest neutral scalar 

# second neutral scalar 

# third neutral scalar 

# lightest pseudoscalar 

# second pseudoscalar 

# charged Higgs 

# neutralino(l) 

# neutralino(2) 

# neutralino(3) 

# neutralino(4) 

# neutralino(5) 

# chargino(l) 

# chargino(2) 

# ~b_l 

# ~t_l 

# ~b_2 

# ~t 2 
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# 3*3 Higgs mixing 
BLOCK HIGMIX 

-1.27437514E-01 # 

-3.16386555E-02 # 

-9.91341856E-01 # 

-2.01793548E-01 # 

-9.77759539E-01 # 

5.71458423E-02 # 

-9.71101974E-01 # 

2.07328914E-01 # 

1.18218769E-01 # 

# 2*2 Pseudoscalar Higgs mixing 
BLOCK AMIX 

1 1 -2.37325986E-02 # 

1 2 -9.99718342E-01 # 

2 1 9.99718342E-01 # 
2 2 -2.37325986E-02 # 

# Gaugino-Higgsino mixing 
BLOCK NMIX # 5*5 Neutralino Mixing Matrix 



1 


2 


1 


3 


2 


1 


2 


2 


2 


3 


3 


1 


3 


2 


3 


3 



.(1,1) 

.(1,2) 
.(1,3) 
.(2,1) 
.(2,2) 
.(2,3) 
.(3,1) 
.(3,2) 
.(3,3) 



.(1,1) 
.(1,2) 
.(2,1) 
.(2,2) 



1,1) 
1,2) 
1,3) 
1,4) 
1,5) 
2,1) 
2,2) 
2,3) 
2,4) 
2,5) 
3,1) 
3,2) 
3,3) 
3,4) 
3,5) 
4,1) 
4,2) 
4,3) 
4,4) 
4,5) 
5,1) 
5,2) 
5,3) 
5,4) 
5,5) 

BLOCK UMIX # Chargino U Mixing Matrix 
1 1 4.31428408E-02 # U_(l,l) 

1 2 -9.99068914E-01 # U_(l,2) 

2 1 9.99068914E-01 # U_(2,l) 
2 2 4.31428408E-02 # U_(2,2) 

BLOCK VMIX # Chargino V Mixing Matrix 



1 1 1.05202671E-01 # N. 

1 2 -7.87589631E-02 # N. 

1 3 6.92565025E-01 # N. 

1 4 -6.93980327E-01 # N. 

1 5 1.46541512E-01 # N. 

2 1 -3.48143940E-02 # N. 
2 2 3.76738427E-02 # N. 
2 3 7.01102969E-01 # N. 
2 4 7.06807991E-01 # N. 

2 5 7.90299289E-02 # N. 

3 1 -2.45218320E-02 # N. 
3 2 9.59570918E-03 # N. 
3 3 -1.58479313E-01 # N. 
3 4 4.52366830E-02 # N. 

3 5 9.85973910E-01 # N. 

4 1 9.93510538E-01 # N. 
4 2 1.73884449E-02 # N. 
4 3 -5.24484443E-02 # N. 
4 4 9.87396900E-02 # N. 

4 5 1.15796113E-02 # N. 

5 1 -7.47306864E-03 # N. 
5 2 9.95983570E-01 # N. 
5 3 3.06884365E-02 # N. 
5 4 -8.37728195E-02 # N. 
5 5 -1.10280322E-03 # N. 



1 1 1.18343118E-01 # V. 

1 2 -9.92972762E-01 # V. 

2 1 9.92972762E-01 # V. 
2 2 1.18343118E-01 # V. 



1,1) 
1,2) 
2,1) 
2,2) 



Table 2: Corresponding spectr.dat output file. 
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# 

# Total number of points scanned 
# 

1000 
# 

# Output format O=short l=long (not recommended for big scannings) 
# 


# 

# lambda 
# 

0.5 
0.5 
# 

# kappa 
# 

-0.15 
-0.15 

# 

# tan (beta) 
# 

3.5 
3.5 
# 

# mu 
# 

200. 
200. 

# 

# A_lainbda 
# 

780. 
780. 
# 

# A_kappa 
# 

150.0 
250.0 
# 

# Remaining soft terms (no scan) 
# 

mQ3= 1.D3 

mU3= 1 . D3 

mD3= 1.D3 

mL3= 1.D3 

mE3= 1 . D3 

AU3= 1 . 5D3 

AD3= 1 . 5D3 

AE3= 1 . 5D3 

mQ= 1.D3 

mU= 1 . D3 

mD= 1 . D3 

mL= 1 . D3 

mE= 1 . D3 

Ml= 5.D2 

M2= 1.D3 

M3= 3.D3 



Table 3: The scaninp.dat file for parameter scan #2. 
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Parameter Scan n 1 
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Figure 1: Branching ratios of h\ as a function of A\ for \ = k = 0.3, tan/3 
/Xefif = 180 GeV, A^ = 0, rngquark = 1 TeV, and At = 1.5 TeV. 
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Figure 2: rrih^ and r/ia^ as a function of Ax for the same parameters as in fig. 1. 
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Parameter Scan n 2 
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Figure 3: Branching ratios of hi as a function of A^ for A = 0.5, n = —0.15, tan/3 
= 3.5, /ieff = 200 GeV, Ax = 780 GeV, m^quark = 1 TeV, and At = 1.5 TeV. 
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Figure 4: Branching ratios of /i2 as a function of A^ for the same parameter choices 
as in fig. 3. 
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Parameter Scan n 2 
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Figure 5: rrih^ and m/^ as a function of A^ for the same parameters as in fig. 3. 
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Figure 6: LEP constraints in comparison to predictions for hi for the parameters of 
fig. 3. Note the correlation of m = rrih-^ with A^ given in fig. 5. 
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Appendix A 

1 General Conventions for the NMSSM 

Below, we define our conventions for tlie tree level Lagrangian of the NMSSM. The 
superpotential for the Higgs fields, the quarks and the leptons of the 3rd generation 
is 



W = htQ- Huf'k - hQ ■ Hd% - KL ■ HdL% + XSH^ -Ha + ^KS^ . (A. 1) 

Hereafter, hatted capital letters denote superfields, and unhatted capital letters the 
corresponding (complex) scalar components. The SU{2) doublets are 



Products of two SU{2) doublets are defined as, e.g.. 



Hu- Hd = H^H^ - H^H^ . (A.3) 

For the soft SuSY breaking terms we take 

--Csoft =mljH^\^ + mljHd\^ + ml\S\^ + mim+ml\Tl\ 
+ml\Bl\ + ml\L'^\ + m2|L|| 
+{htAt Q ■ H^T^ - hAb Q ■ HdB'j^ - KA^ L ■ HdL^^ 

+\AxHu-HdS + \KA^S^ + h.c.). (A.4) 

2 Higgs Sector at Tree Level 

For completeness, we list here the Higgs potential, tree level Higgs masses and our 
conventions for the mixing angles. The tree level Higgs potential is given by 

+Xk{H^ ■ HdS*' + h.c.) + \ /(|/JJ2 - \Hd\^f 

+{\AxH^-HdS+^KA^S'' + \i.c.) (A.5) 
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where 
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9^ = 1 (gf + gl) ■ (A.6) 



Assuming vevs h^, hd and s such that 

n^- liu^ y= , iirf - ftd H y= , b - s^ -= — (A.7j 

eq. (A. 5) simphfies to 

V = X^ihl-s^ + hls^ + hlhl) + /t^/ - 2\KhuhdS^ - 2\Ax KhdS 

+ - kA^s^ + m^>^ + ml hi + m|s2 + _ ^2(^2 _ ^2^2 ^ ^^ g^ 

The sign conventions for the fields can be chosen such that the Yukawa couphngs 
A, ht, hb, the vevs hu, hd (and hence tan/?) as well as the soft gaugino masses M, 
are all positive. Then, the Yukawa coupling k,, the trilinear soft terms Ai, and the 
vev s (and hence fieff) can all be either positive or negative. 

2.1 CP-even neutral states 

In the basis 5'''"'"^ = {Hur, HdR, Sr) and using the minimization equations in order 
to eliminate the soft masses squared, one obtains the following mass-squared matrix 
entries: 

A^lii = g^hl + Xs-^{Ax + Ks), 

-^1,22 = g^hl + \s^{Ax + Ks), 

rid 

h h 
>1|33 = \Ax^^ + KS{A^ + 4/ts), 

-^1,12 = {2\^-g'')Khd-\s{Ax + ns), 

A<|^i3 = 2\^huS- \hd{Ax + 2Ks), 

•^1,23 = 2\^hdS-\hu{Ax + 2Ks). (A.9) 

After diagonalization by an orthogonal matrix Sij one obtains 3 CP-even states 
(ordered in mass) hi = SijSj"'^'^, with masses denoted by m/j.. 

In the MSSM limit (A, /t -^ 0, and parameters such that h^ ~ Sr) the elements 
of the first 2x2 sub-matrix of Sij are related to the MSSM angle a as 



Sii ~ cosa , S21 ~ sina , 

5*12 ~ —sin a, 5*22 ~ cos a. (A. 10) 
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2.2 CP-odd neutral states 

In the basis P^"*"^ = [Hui, Hdi, Sj) and using the minimization equations in order to 
ehminate the soft masses squared, one obtains the following mass-squared matrix 
entries: 

Ml^n = Xsj^{Ax + Ks), 

MI22 = \s-^{Ax + Ks), 
hd 

h h 
A1p33 = AXnhuhd + \Ax— — —SkA^s, 

Mp 12 = As(Aa + /€s), 

•^p,i3 = Xhd{Ax-2Ks), 

M%23 = Xhu{Ax-2Ks). (A.ll) 

The diagonalization of this mass matrix is performed in two steps. First, one rotates 
into a basis {A, G, Si), where G is a massless Goldstone mode: 

COS/3 -sin/3 0\ / A \ 
sin/3 cos/5 G (A.12) 

I J \Sr J 

where tan/5 = hu/hd- Dropping the Goldstone mode, the remaining 2x2 mass 

matrix in the basis {A, Sj) has the matrix elements 

Ml,, = \s^^^^{Ax + ns), 

A1p22 = -iXKhuhd + \Ax— — — SkArS, 

M%,2 = \^hl + hj{Ax-2ns). (A.13) 

It can be diagonalized by an orthogonal 2x2 matrix P/- such that the physical 
CP-odd states Oj (ordered in mass) are 

ai = P[,A + P[2Si 

= P[, (cos PHui + sin pHdi) + P[2Si, 

a2 = P^,A + P^2Si 

= P^,{cos f3H^j + sin f3Hdi) + P^2Si, (A.14) 

and, for completeness, 

G = - sin pH^j + cos pHdi . (A. 15) 
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The decomposition of the bare states in terms of physical states reads 

Hui = cosP{P[iai + ^21*^2) — sinPG , 
Hdi = sin l3{P{^ai + P2i'^2) + cos (3G , 
Si = P[^a, + P^^a2. (A.16) 

(In principle, since the matrix P/ is orthogonal, it could be parameterized by one 
angle.) Eqs. (A.16) suggest the introduction of a 2 x 3 matrix Pij with 

P,i = cos/3/^'i, P,2 = sin/5P^i, P,3 = P:2 (A.17) 

such that, omitting the Goldstone boson, 

Hui = PiiOi + P2ia2 , 

Hdl = -PllOl + -^21^2 , 

Si = Pi3ai + P23a2. (A. 18) 

2.3 Charged states 

In the basis (-f^^, [H^]* = H^l), the charged Higgs mass matrix is given by 

Ml = ^Xs{A, + ^s) + Khdi^j - A^)) ( """"l^ J^^ ) . (A.19) 

This gives one eigenstate H^ of mass TrA^^ and one massless goldstone mode G^ 
with 

H^ = cos(3H^-sm(3G^ , 

H^ = sinpH^ + cos I3G^ . (A.20) 

3 SUSY Particles 
3.1 Neutralinos 

Denoting the f/(l)y gaugino by Ai and the neutral SU{2) gaugino by X^, the mass 
terms in the Lagrangian read 

£ = ^M.X.X^ + ^M^XlXl 

+ ^ Ai(/i„^:! - hd^d) - ^ ^lihui^u - hdi^'d)- (A.21) 
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In the basis ip^ 



-i\i, —i\2, ipuj ipdy i^s) one can rewrite 
1 



C 



i^yMoir) + h.c. 



(A.22) 



where 



Mr 



( Ml 





gihu 

V2 


aihd 

V2 





v 


M2 


g2hu 



92hd 






-\hd 

2ks 



J 



(A.23) 



(Recall that here /i = /icg = As). One obtains 5 eigenstates (ordered in mass) 
Nijip^, with Nij real, with masses m^o that are real, but not necessarily 






Xi 



positive. 



3.2 Charginos 

The charged SU{2) gauginos are A~ = ^ (A2 + iXl), A"*" = 75 (A2 — iXl) 
Defining 



^^ 



the Lagrangian can be written as 



C 



with 



-i^ii'^.r) 



X 



i,- 



x^ 

X 



M2 g2K 
92hd /i 



-iX- 






h.c. 



The mass eigenstates are y^ = Vip'^, x = Uip , with 



U 



cos 9u sin 61 
— sin 9u cos 9i 



V 



cos Oy sin 9\ 
— sin 6v cos 6 



V 



Defining 
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v/Tr(X^X) - 4det(X^X) 



(A.24) 



(A.25) 



(A.26) 



(A.27) 



(A.28) 



one has 



tan 6^(7 
tan6'y 



2g2{M2K + fihd) 
2g2{M2hd + fiK) 



(A.29) 
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where —-njl < 6u, Oy < 7r/2 are such that Md = UXV^ is diagonal, but not 
necessarily positive. The masses with |m^J < jm^jl are given by 

m^^ = cos9u{M2 cos9v + g2huSm9v) + sin 9u{g2hci cos 9v + fj,sm9v), 
m^^ = sm9u{M2 sin^y — g2hu cos9v) — cos 9u{g2hd sin 9v — ficos9v)- 

(A.30) 

/ v+ \ 
In terms of 4 component Dirac spinors ^j = I _*_ 1 one can rewrite the Lagrangian 

as 

£ = -x'MdX^ + h.c. = -m^^,"^!^! - m^^'^2'^2 ■ (A.31) 

3.3 Top and Bottom Squarks 

To complete the consequences of our conventions above, we give here the top and 
bottom squark mass-squared matrices (without the D-term contributions). Below, 
ti, t% bi and 6^ denote the two component quark spinors. 
Top squarks: 





Tr 


Tl 


Tr 


f m\. + h'^hl 


ht{AtK - >^shd 


Tl 


\ ht{AtK - As/id) 


^Q + h^thl 



(A.32) 



Bottom squarks: 



Br Bl 

B*R f m| + hlhl h{Abhd - XsK) \ (A.33) 

Bl I hb{Abhd - Xshu) ml) + /i^/i^ 
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Appendix B 

1 Feynman rules for the Higgs Couplings 
1.1 Higgs- Quarks 

The couplings are obtained by expanding the quark mass matrices in the (properly 
normalized) physical Higgs fields /ij, aj and H^ . Below, we use w^ = /i^ + /i^, and 
consider the quarks of the third generation. 

hf f- ■ "^* c 

yzf smp 



V2v cos (3 



hibLb^ : -j= -Si2 



ailLljl ■ ^ r- . -Til 

V^t'smp 

aibLbn : t^= -Pi2 

y ivcosp 

H^bLfr, : -— cot/3 

V 

H-tLb% : -— tan/5 (B.l) 



t" 



1.2 Higgs- Gauge Bosons 



These couplings are obtained from the kinetic terms in the Lagrangian: 

-\- 
hiZ^Zy : g^y -= — [huSn + hdSi2) 

h^W+W- : g^y^{KSa + h,S,2) 
hi{p)H+{p')W- : ^{cosf3Sa-smpS,2){p-p'), 

a,{p)H+{p')W- : 2^(cos/5P,i + sin/5P,2)(p-p')M 

n 
hi{p)aj{p')Zf, : i^={SiiPji~ Si2Pj2){p-p')f, 

v2 

H^{p)H~{p')Z, : ^/Z^ (p-pO. (B.2) 

1.3 Higgs-Neutralinos/Charginos 

As in the case of the Higgs-Quark couplings, these couplings are obtained by ex- 
panding the corresponding mass matrices: 
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haXtXj 


ttaXtxJ 


H-X7X^ 


H-Xtx^ 


haXh', 






Sa3Ui2Vj2 H -i={SalUiiVj2 + Sa2Ui2Vjij 



I ( —P,^U,^V,2 - ^iPalUaV,2 + Pa2U^2V,^: 



sin 3 

A COS/?f/j2iVj5 ■J=-Ui2{giNjl + fi'2iVj2) + 5-2 sill /Sf/ji A/'j4 



cos/3 
72 

^TT ('S'ain-- — S'a2n- • ) + — (S'aill-- — S'a2n-- 



A sin[3Vi2Nj^ + —^Vi2{giN.j^ + ^2A^i2) + 92 cos f3VnN,^ 



-^{Sal^ij - Sa2^ij ) + —{Sa\^ij - Sa2S^ij ) 



+|(i^ain^| - Pa2n^/) - f (Pain^l - Pa2n?/)) (B.3) 

where ng* = iVi^AT^b + Ni^N^a- 

1.4 Triple Higgs Interactions 

The trihnear Higgs self-couphngs are obtained by expanding the scalar potential. 

^^A-pri23 I ^^« 77333 



a6c ^^ „ /7^ afec 



,2 



where 



+^ (/^«(niL^ - n-) - /.,(n-,^ - n-^) (b.4) 

'rJbii^cjJak + Jci^aj^bk + ^ci^bj^ak ■ 1^.5) 
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haat,a, : ^ {KiU]^l + UlU) + h.iUlU + U^,) + s{Illll + U^J) 



\k 



313\ 






333 



A/TT123 I tt213 I tt312\ '^^k ,-.333 






where 



' /pr \ abc ' abc ' abcJ r^ 

+^ {huiKll - Kfc) - h,{nlll - nllD) (b.6) 



^Zc - Sai{PbjPck + PcjPbk) ■ (B.7) 



+^ {huiK'' - K'') + h,{Uf' - Uf')) (B.8) 



where 

nf = 25,,C,Cfc (B.9) 

with Ci = cos/3, C2 = sin/9. 
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Appendix C 

In this appendix we give the exphcit expressions of the radiative corrections to the 
Higgs masses. 

Generally, the radiative corrections to the Higgs effective action consist in cor- 
rections to the wave function renormalization constants Zj, and corrections to the 
effective potential \4//(^j)- We will treat the one loop corrections originating from 
stop and sbottom loops exactly in the stop and sbottom mixings and mass splittings. 
The resulting contributions to the Higgs masses are nevertheless quite simple if they 
are expressed in terms of couplings and masses at the scale Q = Msusy (the average 
of the squark masses, see below) and in terms of the Higgs vevs before rescaling by 
Z^ , which we will denote by hi{Q). However, the Higgs vevs are related to the 
Fermi coupling Gp and our definition of tan /9 a/ter rescaling by Z^ , cf. eq. (2.9). 

Hence, we first have to determine Z-i for i = H^, iJ^: 

3ft| 

where 

t = In(gV^top) , 
Q' = \{2ml + ml + ml)+ml^ (C.2) 

(such that t = for vanishing SUSY breaking squark masses niQ, rriT and m^). 
Now we have 

K{Q) = KI^ZZ , hd{Q) = hJ^zZ . (C.3) 

with hu and h^ given in terms of Gp and tan 13. 

In addition we need the running Yukawa couplings and quark masses at the scale 
Q. First, the running bottom quark mass (and hence /i?,) at the scale mtop is given 
by its QCD evolution to one loop order 

hbirutop) = hirribot) 1 - -^^-"^("^top) In —^ . (C.4) 

The running top quark Yukawa coupling at the scale mtop is given in terms of the 
top quark pole mass as " 

mt op{pole) f Aasimtop) llalim top) 
hu V Stt vr^ 



ht{mtop) = "^^T " 1 + ^^^^F^ + ^ . ^^ • (C.5) 



"Here we use the MS relation between the pole mass and the running mass. In the DR scheme 
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Next, the running Yukawa couplings at the scale Q are obtained from the RG 
equations (for an effective 2 Higgs doublet model) 



i9ht + hi) - ht^ , 



dlnQ2 647r2' ' "' n 

dhb hb 2 t2^ ■ <^ 



i9hi + hi) - h— . (C.6) 



dlnQ2 647r2' " '' n 

Its solutions, exact in as but perturbative in the Yukawa couplings, are 

ht{Q) = htimtop) (l + J;^as{mtop)t\ (l + -^{9h^t+hl)t 

h{Q) = h{mtop) (l + J;^as{mtop)t\ (l + -^{9hl + h^^)t\ . (C.7) 

However, possibly a SU{2) Higgs doublet has a mass much larger than rritop- 
Then it is practically degenerate, and its approximate masses are given in the 
NMSSM by 

MHD = ^s{Ax + Ks){coil3 + iaiil3) . (C.8) 

At scales below Mhd the RG equations read 

dh^ = ^(9^*+3M-/^.-, 
dh hb ,^,2 , „,2^ • " 



In order to correct the expressions for ht{Q), hb{Q) in this case, we have to mul- 
tiply them by (^1 + ^\n{M]jjj/mj^p)j , (^1 + ^\n{Mfjj^/mf^p)j , respectively, if 
Mhd > rritop- 

The expressions for the running top and bottom quark masses at the scale Q are 
then simply 

mtopiQ) = K{Q)K{Q) , mbot{Q) = hb{Q)ha{Q) , (C.IO) 

and below it is convenient to use 

tan(3{Q) = K{Q)/hd{Q) . (C.U) 

For the (s)top, (s)bottom induced one loop corrections to Vg// we still need the 

eigenvalues of the stop, sbottom mass matrices (A. 32) and (A. 33). To this end it is 

the second factor 4/37r would read S/Stt, which would increase ht{mtop) and hence the lightest 
Higgs mass by ^ 1%. Since we have not included subdominant (single) logarithms in the two loop 
corrections to the effective potential we are not sensitive to the scheme in which the running top 
quark mass is defined, which leads to a theoretical error of ~ 1% on the mass of the lightest Higgs. 
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convenient to define, step by step, 



Mlop = 2K + ^t) 



AMlp = -{ml -ml) , 

AMI, = lirnl-ml) , 

Xt = At- \s cot P{Q) 

Xh = Ah- Xstd.n(3{Q) 



Wt = JAM^top + rnl^{Q)X! , 



^b = \I^M%^, + mUQ)Xl, (C.12) 

with the help of which the eigenvalues of the stop, sbottom mass matrices become 





= Ml^^ + ml^{Q)-Wt, 


^^^stop2 - 


- Ml^^ + ml^{Q) + Wt, 




- Mi, + mL(g)-w-6, 


Mlot2 -- 


= Ml^^ + ml^{Q) + Wh. 



(C.13) 



Finally the following auxiliary quantities are useful: 



'"'top^t / ^^^stop2 ^ ^^^stopl , / ^^^stop2 



„jy.f — "-"y ■- I ai-cy^ ai-cpi 1 1 !,Lupz. 1 _ 2 

Wt' \M^top2-MLpi \MLpi 

, ^top^t ( ^stop2 

emt = — -t: — In ' 



Wt "'\M^topi 

nmh - ^L^6 f MLt2 + Ml,, /Ml,, ^ 

^^b \^^^sbot2 ^^^sbotl \^^^sbotl 



mLXh, /M|, 
r2 



^"^^"rnA^i' <^'"' 



and, remarkably, quite many of the radiative corrections can be absorbed into a shift 
oiAx: 



Subsequently the index new of Ax^new will be omitted. 



AA,„e«, = Aa + -^^At fmt + ;^7 Afc /m6 . (C.15) 
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Then the radiatively corrected mass-squared matrix entries in the CP even sector 
are, instead of (A.9), 



-Mill = g^K{Qf + \s{A^ + Ks) cot P{Q)-^-^^^^{\sf gmh 

+ oo 2 ^^t gmt + AAt emt + Am^^^ In ^ 

6Zn \ \ ITT'top 

-M522 = g'haiQY + Xs{A^ + Ks) ta.n (3{Q) - ^^{XsY gmt 



.!^ (-^f <,,„, + 4^, em6 + 4mL In ( "'^^f^" 



^bot 



-^^A^/^.(g)^ ^mt - ^^A^/.„(g)^ gmb , 
-Ml 12 = (2A2-/)/^„(Q)/i,(Q)-As(A;, + «:s) 

A<|i3 = 2A2/i„(g)s-A/id(g)(AA + 2fi:s 



+As ( ^*^_2 (^* ^^^ ^ 2 emt) + ^t;_; (^b gmb - 2 em6) 



+^^^A/irf(g)(A gmt - 2 emt) 



f ,!;_2 ^'^ ^«(^)(4 M^ - 9mh) , 



327r 

A<|23 = 2X^hd{Q)s - Xhu{Q){Ax + 2Ks 
, 3/i,(g)2 



327r2 
327r2 



A/;.M(g)(A6 f?^ — 2 em6) 

A^s /id(g)(4 /mt - gmt) . (C.16) 



The dominant two loop effects ~ h^t"^ and ~ hjast^ contribute to A^ln only. They 
read 

5Ml,, ,,^^^, = ^-^^KiQr (327ra,(g) - hm'^ t' . (C.f7) 

The above expressions for the elements of the CP even mass matrix have been 
obtained by taking the second derivatives of the effective potential, and using the 
three minimization equations with respect to /i„, h^ and s in order to eliminate the 
three SUSY breaking masses squared. Since the Higgs fields are not quite properly 
normalized due to the Z factors, the above mass matrix elements still have to be 
rescaled by appropriate powers of Zi (here below Ai'sij norm denote the mass matrix 
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elements for correctly normalized Higgs fields): 



2 

5,11 norm. 

2 

5,22 norm. 




2 

5,12 norm. 


= ^^u/V^H^^H^ 


2 

5,13 norm. 


= MI,J^Zh^ , 


2 

5,23 norm. 


= MI,J^Zh,. 



(C.f8) 

Now we turn to the electroweak leading logarithmic corrections to the Higgs 
mass matrix elements, which are taken from ref. [15]. As in [15] we allow for a 
Higgs doublet which is possibly heavy {Mhd > Mz), and for gauginos/Higgsinos 
with masses possibly (much) larger than Mz- In this case their masses are given 
approximately by the soft SUSY breaking gaugino mass parameters Mi, M2 or 
Heff = As, respectively. The following auxiliary parameters are useful subsequently: 

Hi = max{\fieff\, Ml} , 
H^ = max{\fXeff\, M2} , 
/i^2 = max{\i2eff\, Ml, M2} . (C.19) 

The electroweak corrections to the Higgs mass matrix elements are then given 
by (we omit the index norm, in the following) 
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(sin^ Ow + cos^ Ow) Hfilff/Ml] 



IQlX^My^ \ 3 

+ -[-9(cos2 13 - W (3Y + (1 - 2 sin^ ^vk + 2 sin^ dw){coi^^ P - sin^ P? 
6 

+10 - 2 sin^ ew + 2 sin^ ^vk] ln(M^o/M|) 

+ sin^ ewil-2 sin^ ^vf) ln(/i?/M|) - 4 sin^ ^vf cos^ Ow Ivl^^iI^/MI) 

+ cos^ Owi'i - 10 cos^ ^vf) ln(/i^/M|) - - cos^ ^vk ln(M|/M|) 

o 

2 

-3' 

5e»A^|,22 = Y^5^cos='/5((4-8sin2^p^ + — sin^^H^)t 

+ -[-9(cos2 /3 - W /3)^ + (1 - 2 sin^ ^h- + 2 sin^ ^vk)(cos2 P - sin^ /3)^ 
o 

+10 - 2 sin^ ^p^ + 2 sin^ Ow] ln(M|^/M|) 

+ sin^ ewil-2 sin^ ^vf) ln(/ii/M|) - 4 sin^ ^vf cos^ ^vf ln(/ii2/M|) 

+ cos^ ^vk(3 - 10 cos^ 9w) Hfil/Ml) - ^ cos^ 9w ln(M|/M|) 

Q^M^ ( 32 

5e«;A^|,i2 = Y|^5^sin/5cos/?( (-4 + 8sin2^p^-ysin^^iy)t 

-i[-9(cos2 /3 - sin^ pf + (1 - 2 sin^ ^^^ + 2 sin^ ^vf)(cos2 /? - sin^ pf 
6 

+8 - 22 sin^ Qy^ + 10 sin^ Qw] \n{MJj^/Ml) 

- sin^ ^^(1 + 2 sin^ Qw) ln(/i?/M|) - 4 sin^ Qw cos^ ^p^ ln(/i^2/M|) 

4 

- cos^ ^vf(3 + 2 cos^ ^vf) ln(/i^/M|) + - cos^ Bw HM^/M^) 

+ '^{sm^9w + cos^ew)Hf^lff/Ml)] . (C.20) 

The eigenvalues of this Higgs mass matrix give the "running" Higgs masses at 
a scale rritop due to the argument of the logarithm t in the Z factors (C.l). In the 
LLA we can identify the scales rritop, Mz (in the arguments of the logarithms) and, 
to start with, M Higgs-, and take these "running" masses as pole masses. However, 
if a Higgs mass is (much) larger than m^op, we should have computed the Z factors 
(C.l) at a scale Mfjiggs rather than rritop, in order to identify the eigenvalues of 
the Higgs mass matrix with the Higgs pole masses. We can correct for this at the 
end (after having diagonalized the Higgs mass matrix), taking care that now the 
eigenstates hi are given by hi = SijS^^^, cf. the notation below eq. (A. 9). Hence, if 
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Mmggs > iTT-top, the Higgs pole masses Mf{pole) are given in terms of the eigenvalues 
Mf as 

M^ipole) = Mfl [l - {-^^Sl + ^5|) In {Ml/ml^)^ . (C.21) 

Now we turn to the radiative corrections to the CP odd mass matrix. After the 
shift (C.15) in Ax these are trivial: It suffices to take the tree level expressions (A. 11) 
for the mass matrix elements, to replace the Higgs vevs h^ and h^ by its expressions 
(C.3) at the scale Q, and to rescale the mass matrix elements by appropriate powers 
of Zi as in (C.18) for the CP even sector. After dropping the Goldstone mode this 
gives a 2 X 2 mass matrix (instead of (A. 13)) 

Mp^n = M^x + t^s) — + —— , 

V ^Ha ^Hu J 

Ml,, = 4A«:/.„(g)Mg) + AA.M^IM^ - 3/.A.. , 



-^p,i2 = ^\/hu{Qy/ZH, + ha{Qy/ZH^ {Ax - 2ks) . (C.22) 

Finally, if a CP odd Higgs mass is (much) larger than rritop, its pole mass is given 
by an expression analogous to (C.21) after replacing Sn by P/^cos^/9 and Si, by 
i^'i sin^/3 (see (A. 14) for the definition of Plj). 

In the expression for radiatively corrected charged Higgs mass we neglect stop/- 
sbottom loops (without large logarithms), but we include large logarithms propor- 
tional to electroweak gauge couplings: 

Ml, = (As(A, + «:s) + MQMQ)(^2V2-A^))f^^^ + ^^^') 



+?>hlhlt/{lQV2Ti'^GF) 

+ ^^ (l2t + 3tan2^^(ln(/x?/M|) +41n(/iyM|)) 

-3 H^l/Ml) - 4 HMl/Ml) - 2 \n{^%f/Ml)) (C.23) 



This time, once Mh± is larger than mtop-, its pole mass is given by 
M^4po/e) = M^,/(l-(^cos^/5 + ^sin2/5)ln(M^,/m?j) . (C.24) 
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